Single-photon scattering on a strongly dressed atom 
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We develop the generahzed rotating-wave approximation (GRWA) approach (Phys. Rev. Lett. 
99, 173601 (2007)) to study the single-photon scattering on a two-level system (TLS) with arbi- 
trarily strong coupling to a local mode in a one-dimensional (ID) coupled resonator array. We 
calculate the scattering amplitudes analytically, and the independent numerical results show that 
our approach works well in a broad parameter region. Especially, when the resonator mode is far- 
off resonance with the TLS, we obtain more reasonable results than the ones from the standard 
adiabatic approximation. Our approach is further extended to the cases with a ID resonator array 
strongly coupled to more than one TLSs. 
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I. INTRODUCTION 

Recently, many attentions have been paid to the pho- 
ton transport in a low-dimensional array of coupled res- 
onators |l|-[25| or a one-dimensional (ID) wave guide [24- 
l49j . which plays crucial roles in the realization of all- 
optical quantum devices. In these systems, the two-level 
or multi-level devices coupled to the resonators or wave 
guides can be used as quantum switches to control the 
scattering or transport of the photons. 

Up to now, the single-photon scattering amplitudes in 
the ID resonator arrays or wave guides coupled to a sin- 
gle two-level system (TLS) EHl, [H, [28l - [33j or a single 
three-level s yste m p^ . l4ll . [47ll48| or multi quantum de- 
vices [§-[1^ IS^ - lssf have been well- investigated by dif- 
ferent authors. The relevant multi-photon scattering am- 
plitudes [2l| - [27l m-Hil are also studied. To our knowl- 
edge, all these studies are based on the rotating-wave 
approximation which is applicable under the conditions 
that the frequencies of the two-level or multi-level system 
are very close to the photon frequencies in the resonators 
or wave guides, and the intensities of the coupling be- 
tween the two-level or multi-level system and the photons 
are much smaller than their frequencies, namely, in the 
weak-coupling limit. 

In the experiments, the ID resonator arrays or wave 
guides can be realized with photonic crysta ls ISOl ISll . 
superconducting transmission line resonators |12Hl4l |29| . 
1471 [5^ or other solid devices, while the two-level or 
multi-level systems can be implemented with either nat- 
ural atoms or the solid-state artificial atoms. In the hy- 
brid systems of solid-state devices, it has been predicted 
that -56] one can realize ultra-strong TLS-photon cou- 
pling with intensities comparable or even higher than the 
photon frequencies. On the other hand, the frequency 
of the solid-state TLS can also be controlled easily in a 
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broad region. Therefore, in the quantum network based 
on solid-state devices it is possible to reach the strong- 
coupling and far-off resonance regime where the condi- 
tions of rotating-wave approximation are violated. In 
these parameter regions, we need to develop new meth- 
ods for the scattering between the flying photon and the 
TLS, and then investigate the possible new effects given 
by the strong TLS-photon coupling to the photonic trans- 
port. 

In this paper, we study the single photon scattering on 
the TLS, which strongly couples to the local mode of the 
resonators. Beyond the rotating-wave approximation, we 
develop an analytical approach for the approximate cal- 
culation of the single-photon scattering amplitudes in a 
ID single-mode-resonator array coupled to a single TLS. 
Our method is based on the generalized rotating-wave ap- 
proximation (GRWA) [57] developed by Irish for the sys- 
tem with a single TLS coupled to a single-mode bosonic 
field. 

With our GRWA approach we obtain the single-photon 
scattering amplitudes under the condition that the pho- 
ton hopping between different resonators is weak enough. 
We will show that our approach works significantly well 
in a very broad parameter region, including the region 
where the rotating-wave approximation is applicable and 
the one where the TLS-photon coupling is strong while 
the frequency of the TLS is close to or smaller than the 
photon frequency. 

Especially, in the far-off resonance region where the 
TLS frequency is much smaller than the flying photon, 
the standard adiabatic approximation does not work well 
in the current hybrid system, while our approach still pro- 
vides good results. On the other hand, when photon- TLS 
coupling is strong enough, our GRWA approach shows 
that the single-photon scattering by the TLS becomes 
equivalent to the transport of a single photon in a ID 
resonator array in which the frequency of a certain res- 
onator is shifted. Then our results are furthermore sim- 
plified and one can make reasonable qualitative estima- 
tions for the photon transport characters, even without 



the quantitative calculations. We also show that, our 
GRWA approach can be generalized to the systems with 
a ID single-mode-resonator array coupled with more than 
one TLS. 

This paper is organized as follows: In Sec. we de- 
velop the GRWA approach for the hybrid system of a 
ID resonator array coupled to a single TLS. In Sec. IIIIl 
we analytically calculate the single-photon scattering am- 
plitudes in such a system with our GRWA approach 
and compare our results with the numerical results. In 
Sec. HVl we discuss the single-photon scattering problem 
in the cases of strong TLS-photon coupling. In Sec.|Vl we 
show the GRWA approach in the system with more than 
one TLS. There are several discussions and conclusions 
in Sec. |VD 
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FIG. 1: (Color online) Schematic configuration for the hybrid 
system of a ID resonator array interacting with a TLS. The 
frequency of the photon in each resonator is u, while the in- 
tensity of the photonic hopping is ^. The TLS with frequency 
f2 is localized at the 0-th resonator and coupled to the photon 
with coupling strength A. 



II. THE GRWA FOR THE TLS-COUPLED ID 
RESONATOR ARRAY 

In this paper we consider the transmission of a single 
photon in a ID single-mode-resonator array coupled to a 
TLS which is located in a specific resonator. To obtain 
the reasonable analytical results beyond the rotating- 
wave approximation, in this section we will generalize 
the GRWA approach proposed by Irish in ReL [53| for a 
single resonator coupled to a TLS to our current hybrid 
system of resonator array. We will first show the adia- 
batic approximation in our system for the cases where 
the frequency of the TLS is much smaller than the pho- 
ton frequency, and then develop the GRWA approach as 
an improvement of the adiabatic approximation. For the 
reader's convenience, in appendix A we introduce Irish's 
GRWA approach for the system with a single resonator 
in the view of adiabatic approximation. 



is satisfied. Therefore the term ar^a'-j^^ + h.c. has been 
neglected in our consideration. This condition also pro- 
vides us a small parameter ^/w which is very useful in 
the following calculation. 

The Hamiltonian Ha of the TLS and the interaction 
Hj between the TLS and the photons in the 0-th res- 
onator are 

Ha - ^a. (3) 

and 

Hi = Xa^{al + ao) (4) 

respectively. Here fl is the energy spacing between the 
ground state \g) and the excited state |e) of the TLS, A is 
the coupling intensity and the Pauli operators CTz and ax 
are defined as CTz — |e) (e| — |g) {g\ anda^ = |e) {g\ + \g) (e|. 



A. The system and Hamiltonian 

As shown in Fig. [U we consider a ID array of infinite 
identical single-mode resonators with a TLS located in- 
side a certain resonator, which is marked as the 0-th res- 
onator in the array. We further assume that the photons 
can hop between neighbor resonators. Then the total sys- 
tem is modeled by the Hamiltonian H = He + Ha + Hj , 
where the tight-binding Hamiltonian He of the resonator 
array is 

+ C30 +00 

Hc=0J ^ a]aj-C ^ {a]aj+i + h.c). (1) 

j — — oo j— — oo 

Here w is the frequency of the photons in the resonators, 
^ is the hopping intensity or the inter-resonator coupling 
strength, aj and a| are the annihilation and creation op- 
erators of the photons in the j-th resonator respectively. 
In this paper we assume the weak-hopping condition 

(2) 



B. The adiabatic approximation for the 
TLS-coupled ID resonator array 

In this and the next subsections, we develop the GRWA 
approach for the TLS-coupled ID resonator array. To 
this end, in this section we first formally develop the adi- 
abatic approximation for our system in the far-off reso- 
nance region with << oj. We will show that, due to 
the energy band structure of the photons, the adiabatic 
approximation cannot be used in our hybrid system, even 
in the far-off resonance case. However, as shown in the 
next subsection, we can also develop the GRWA approach 
as an improvement of this formal adiabatic approach, 
and the intrinsic problem of the adiabatic approximation 
in our system is naturally overcome by our GRWA ap- 
proach. Finally, our approach is applicable in a broad 
parameter region, including the cases of £7 << uj and 

Under the far-off resonance condition il << uj, the 
TLS is considered to be the slowly-varying part of the 
system, while the ID resonator-array is the fast-varying 



one. Then we decompose the total Hamiltonian H as 
H — Hi + H2 where Hi = He + Hj is the Hamiltonian 
of the fast-varying part together with the interaction and 
H2 = Ha is the Hamiltonian of the slow- varying TLS. 

The straightforward calculation shows that Hi has the 
eigen-states 
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:A(fc)t"W|0)® |±). (5) 



Here the TLS states |±) are the eigen-states of the op- 
erator ax with eigen- values ±1, and |0) is the vacuum 
state of the resonator array. The photon momentum k 
can take any value in the region (— 7r,7r], the creation 
operator A{ky for a photon with momentum k is given 

by 
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and 



(6) 
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is the set of all the numbers n (k) . 

The operator U in Eq. ([5]) is determined by 



U~^HiU 
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with C an unimportant c-number. This equation leads to 
the result 
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Loi = {lu + ^/u^ - 4e)/2. 
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The results in Eqs. ()10|lip are derived in appendix B. It is 
pointed out that the parameter aj exponentially decays 
to zero in the limit \j\ — > 00. 

Obviously, the eigen-energy of Hi with respect to the 
eigen-state |±,n) is 

E{±,ri) ^ N{n)uj + 2^^n{k)cosk+C (12) 



with 



Nin) = J2n{k). (13) 
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FIG. 2: (Color online) (a) The energy spectrum of the Hamil- 
tonian Hi of the ID resonator-array and the interaction be- 
tween the resonator-array and the TLS. The low-excited spec- 
trum has clear band structure. The intra-band transitions 
are missed unreasonably in the adiabatic approximation, (b) 
The energy spectrum of the Hamiltonian Hi of a single-mode 
bosonic field and the interaction between the single-mode 
bosonic field and a TLS. In such a simple system there is no 
band structure in the spectrum of Hi, and then the adiabatic 
approximation is applicable when Q, « lj. 



As shown in Fig. [2{a), in the weak-tunneling case with 
^ << uj the low-excited spectrum E{±,n) of the eigen- 
energies of Hi has a clear band structure. Each energy 
band includes all the energy levels with the same total 
photon number N(ji) and different photon-momentum 
distribution n. The A'^-th energy band is centered at Nu 
with band width 4iV^. Therefore, in the low excitation 
cases the inter-band energy gap has the same order with 

Ul. 

The spirit of the adiabatic approximation fos'-'GS'] is 
that, during the quantum evolution, the motion of the 
slowly-varying part of the system would follow the mo- 
tion of the fast-varying one. Then the quantum state 
of the fast- varying ID resonator array would be frozen 
in each adiabatic branch with fixed quantum number n. 
Mathematically speaking, in the adiabatic approxima- 
tion for our system, all the i/2-induced quantum transi- 
tions between the states |±,n) and |±,n') with n ^ n' 
are neglected. This treatment leads to the approximate 
eigen-states of the Hamiltonian H as 



I* 



-,n 



(14) 



Now we point out that, the adiabatic approximation is 
not a reasonable approximation for our present system, 
even in the case offl « ui. That is because, the energy 
spectrum of Hi has a band structure, each band includes 
all the states \i,n) with the same total photon number 
N{n). On the other hand, in the adiabatic approxima- 
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tion all the ff2-induced transitions between the eigen- 
states of Hi with different quantum number n are ne- 
glected. It means that, all the inter-band transitions and 
the intra-band transitions between the states \+,n) and 
\ — ,n') with n ^ n' are omitted. In the case of Q << w, 
the omission of intcr-band transitions are reasonable be- 
cause the energy gap between different bands is of the 
order of w, which is much larger than the intensity fl 
of i?2- However, the neglecting of the intra-band transi- 
tions are unreasonable because the energy gaps between 
the levels in the same band can be arbitrary small. 

As a comparison, we re-consider the simple system 
with a single- mode bosonic field interacting with a TLS. 
In Fig. 2(b), we plot the spectrum of the Hamiltonian 



(15) 



with (a, at) the annihilation and creation operators of the 
bosonic field. It is apparent that in such a simple case 
the spectrum of ffRabii does not have any band structure 
(shown in appendix A). Therefore the adiabatic approx- 
imation 58] is applicable when << w. 

Although the adiabatic approximation cannot be used 
in our system, as shown in the next subsection, we can 
still develop the GRWA approach for our system with the 
help of the basis {|^'±.n)}- Furthermore, in our GRWA 
approach the above intrinsic problem of the adiabatic 
approximation is naturally overcome. 



C. The GRWA approach for the TLS-coupled ID 
resonator array 

Now we develop the GRWA approach for the system 
of the TLS-coupled ID resonator array. Similar to the 
original GRWA for the single-mode bosonic field coupled 
to a TLS (appendix A) , our general GRWA approach can 
also be considered as an improvement of the adiabatic 
approximation. 

In our GRWA approach the total Hamiltonian H is 
approximated as Hq which is defined as 



— X/ X/ ^n'%\'^a,n){'^a,n\SNin),Nin') 
n.n' Q— 

+ ^H'^C-\^ + ,n){'i'-.n'\SN(n}^N{n')~l+h.C., 
n,fi' 

(16) 



where the matrix elements 11^;°^^ are defined as 



(17) 



while the symbol 5a,b is defined as 5a,b = 1 for a = h and 
5a,b = for a 7^ 6. 

Apparently, in the approximate Hamiltonian Hq we 
only take into account the quantum transitions between 
the states \^a,n) and p.n') with 



N{n) -I- iV„ = N{n') + Np 



with Na = 1 for a = + and = for a = — . Obvi- 
ously, this treatment is a direct generalization of the one 
in the development of the GRWA for the TLS-coupled 
single- mode bosonic field (shown in appendix A). 
Furthermore, as proved in appendix C, we have 



for N{n) = N{n). 



(19) 



Then the intra-band transition only occurs between the 
states l^a.ii) and \^a,n') with N{n) = N{n'). Therefore, 
all the intra-band transitions, which are unreasonably ne- 
glected in the adiabatic approximation, are included in 
our GRWA approach through the first term of the right 
hand side of Eq. p^ . Thus the intrinsic problem of the 
adiabatic approximation is overcome in our GRWAap- 
proach. Therefore, our approach is applicable in the far- 
off resonance cases with << uj. 

On the other hand, it is easy to prove that, under 
the near-resonance condition Q ^ uj and weak-coupling 
condition |A| << uj, our GRWAapproach returns to the 
rotating-wave approximation. Therefore, as shown in 
Sec. Ill, our approach is applicable in a broad param- 
eter region with < uj. 

Finally, we point out that, the approximate Hamilto- 
nian Hq in the GRWA approach in our case can also be 
re-expressed in a simple form 

Hg = UHf^^U-^ (20) 

with U given by Eq. © and defined as 



H 



RWA 



[Pn {U- 



'HU) Pn 



(21) 



Here the operator P„ is defined as 

' |0)(0|® |5)(<?|, n = 

mi ,n~ 



Pu = < 



J2{m,}Ui \m)i{mi\ (g) \e){e\Sj2 



^ mi ,n— 1 1 



where \mi)i is the Fock state of the l-th resonator with m/ 
photons. It is apparent that Pn is the projection operator 
to the eigen-space of the total excitation operator 



+ 00 



|e)(e|+ J2 



(22) 



j = -oo 



(18) 



with respect to eigen-value n. In this sense, similarly 
to the GRWA for the TLS-coupled single-mode bosonic 
field discussed in appendix A, our GRWA approach can 
also be considered as "the rotating-wave approximation 
for the rotated Hamiltonian U~^HU". 
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III. THE SINGLE-PHOTON SCATTERING 
AMPLITUDES 

In the above section we generalize the GRWA to the 
hybrid system with a ID resonator array coupled to a 
single TLS. The single-photon scattering amplitudes in 
such a system have been calculated analytically under the 
rotating- wave approximation [2,]. In this section we cal- 
culate the single-photon scattering amplitudes with our 
GRWA approach which is applicable in a broader param- 
eter region. 



A. The single-photon scattering amplitudes 

The single-photon scattering amplitudes can be ex- 
tracted from the asymptotic behavior of the eigen-state 
of the Hamiltonian i7, which is approximated as Hq in 
our GRWA approach. To this end, we need to solve the 
eigen-equation 



HG\^{k))=E(k)\^{k)) 
with boundary conditions 



(23) 



l*(fc)> = (e''=(-^")+rfce-'^(-^))|l)_,|0),|<I>'(fc)) 

+t,e^^^|0)_,|l),|$'(fc)) + |0)_,|0),|$(fc)) 

(24) 

in the limit of j — > +oo. Here |0)_j|0)j is the vac- 
uum state of the resonator modes in the j-th and 
—j-th resonator, |l)_j|0)j and |0)_j|l)j are defined as 
a-j\0)-j\0)j and a]|0)_j|0)j respectively. |$(fc)) and 
|$'(fc)) are the quantum states of the TLS and other 
resonators except the ±j-th ones, and tk are the 
single-photon reflection and transmission amplitudes, or 
the single-photon scattering amplitudes. 

The physical meaning of the boundary condition (j24p 
can be understood as follows: For the scattering state 
with respect to a single photon input from the left of 
the ID resonator array, there are three possible rele- 
vant states for the —j-th and j-th resonators with large 
i.e., \0)-j\0)j, and |0)_j|l)j. Furthermore, 

the possibility amplitude with respect to |l)_j|0)j is 
^g«fe(-i) _l_ ^^g-«fc(-j)^ ^ since the photon in the —j-th 
resonator can be either the input one or the reflected 
one. Similarly, the possibility amplitude with respect 
to |0)_j|l)j is ifce*'^-'. It is easy to prove that the 
boundary condition used in the calculation of the single- 
photon scattering state with rotating-wave approxima- 
tion (Eq. (5) of Ref. 0) can be re-expressed as the one 
in Eq. 

Usually the expression of Hq in Eq. (^0)) is compli- 
cated and the eigen-equation ((23l) is difficult to be solved 
directly. However, due to Eq. I^U^, the Hamiltonian Hq 
is related to H^^^ through a unitary transformation. 
Then the eigen-equation (1^51) of Hq is equivalent to the 



one of 7J|^WA. 

Hr^\^Rik))^Eik)\^nik)) (25) 

and the eigen-state I^^r (fc)) of H]^^ is given by 

\^nik))^U-'\^{k)). (26) 

More importantly, with the help of Eq.Q and the fact 
that lim|j|_^oo aj = 0, the boundary condition ((M)) for 
1 4' (fc)) is transformed to the one of (fc)), i.e., in the 
limit of j oo we have 

\^n{k)) = (e^'^^-^) +r,e-''=(-^)) |l)_,|0),|$i, (fc)) 
+i.e^^-^|0)_,|l),|<i>i,,(fc)) + |0)_,|0),|$fl(fc)), (27) 
where (fc)) and |$'^ (fc)) are defined as 

- Y[ciip[~a,a,ial-a,)]\<i>ik)), (28) 

and 

l^Rik)) = l[cxp[-~a,a,{al-aM^'{k)). (29) 

respectively. Therefore, the scattering amplitudes and 
tk can be obtained from the solution of the eigen-equation 
of H^^^ with boundary conditions (^7)) . 



B. The perturbative approach for the single-phton 
scattering amplitudes 



^RWA 



In this subsection we solve the eigen-equation of 
and calculate the single-photon scattering ampli- 
tudes. To this end, we first use the explicit result about 
the unitary operator U shown in appendix B to calculate 
the Hamiltonian U^^HU. We have the results as: 



U-^HU 



' cosh(^^ 2vi)az — i sinh(^^ 2vi)ay 



(30) 

with Vi = — ai(aj — ai). In principle, we can derive the 
exphcit expression of with Eqs. ((HI) and ((30l) . 

Here, for simplicity, we expand H^^^ as a power series 
of the parameter ^/uj and only keep the low-order terms 
under the weak- hopping condition ^ << uj. Then we can 
analytically solve Eq. ((25)) with the approximated 
and derive the single-photon scattering amplitudes. 

Now we calculate the single-photon scattering ampli- 
tudes with first-order approximation where only the 0-th 
order and 1-st order terms of ^/w are kept in H^^. As 
shown in the following, in most of the cases, this approxi- 
mation is enough to give accurate results for the scatter- 
ing amplitudes. The straightforward calculation shows 
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that, up to the first order of ^/w, is approximated 

as 

rrRWA ^ rrRWA(l) 
R H 

j j 
+4°;(|0<?)(05| +^ +4e'|0e) (Oe| 

+J°J\hg){hg\ + J<^"H\Oe){log\+h.c.) 

+uj['J{\lig){lo9\ + \l-ig){lo9\+h.c.) 

+J« (|0e) (li.g| + |0e) (l_i.g| + /i.e.) (31) 

with the parameters 




In Eq. ([5T|) the states |0e), IO5), |lje) and are 
defined as |0)|e), |0)|(7), |lj)|e) and respectively, 
where |0) is the vacuum state of all the resonators. 
The physical meaning of Eg. (1311) is very clear. In the 0- 

th order terms of £,/uj, or the terms proportional to i^og^, 

WQg\ uj^^J and J^^\ the effective couplings occur between 

I 



the TLS and the photon in the 0-th resonator in which 
the TLS is located. Nevertheless, in the 1-st order terms 
proportional to J^^\ effective couplings appear between 
the TLS and the modes in the ±l-st resonators. These 
terms imply that the non-rotating-wave effects from the 
coupling between the TLS and the 0-th resonator can in- 
directly influence the behavior of the modes in the ±l-st 
resonators. Furthermore, the hopping intensities between 
the 0-th and ±l-st resonators are also tuned by the terms 
with uj[]^ . 

As shown above, the single-photon scattering ampli- 
tudes are approximately derived from the eigen-equation 

Hf^^'^\^nik))=E{k)\^j^ik)) (33) 

of jj^^^^") with boundary condition pT]). It is apparent 
that the solution (fc)) of Ea. (P5|) takes the form 

+00 

\^R{k))^ ^ ui'\j)\l,g}+uW\Oe} (34) 
with the coefficients u), {j) given by 

u^k\j) = \ MO), J=0 • (35) 

Substituting Eqs. (l34l35p into Eq. ([55)) . we obtain the lin- 
ear equations for the reflection amplitude r^}'^ and trans- 
mission amplitude t^i}\ These equations can be solved 
analytically. Then we obtain the scattering amplitudes 
rl^'' and tL^'' given by the first-order approximation of 

^RWA. 



4A2^rj2 cos(A:) -t- 2e^uj^ [w^ -|- 2^uj cos(fc) - 4^^ „ 4^2 (2fc)] 



„(i) 



A^uj (2A2 + cos (k) + {AX^ 



■ uj cos 



(2fc)] 



_4g2»fe;^4^^3 ^ _^ g2»fe-) ^^6 _ 3^ COS {k)] 

+e'^uj^n [8e3''=A2^2 ^ ^e^'^X^ [A^ - UJ^) + (2A2 + uj^) - e^^'^^uj (QX^ + w^)] 



(36) 



2ik+- 



' A2f72 [4ecJ (2A2 + + (cj4 _ 4^2^2 _ IQX2^2-^ _ ^^4 gjj^ 



1. 



(37) 



The above procedure can be straightforwardly gener- 
alized to the cases with high-order approximations of 
H^^. For instance, in the second-order approximation, 
is approximated as ff^^^(^) which includes the 0- 
th order, 1-st order and 2-nd order terms of ^/cj. It can 
be found that in ^^■^^f^) ^-^q tLS is effectively coupled 
to the 0-th, ±l-st and ±2-nd resonators. We can also 



solve the eigen-equation of Hj^ , and obtain the an- 
alytical expressions of the relevant scattering amplitudes 

(2) (2) 

Tj, and t\ . In general, for any integer n, the scattering 

amplitudes r^"'' and t^^^ from the n-th order approxima- 
tion of can be obtained with the similar approach. 

In the limit of rt — >■ cxd, the results r^"^ and t^"' would 
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FIG. 3: (Color online) The single-photon reflection rate \rk\'^ given by our GRWA approach with the first and second order 
approximation for i.e., 1"^ (red solid line) and |r^^' p (black empty triangle), |rfep from the rotating-wave approximation 

(RWA) (green empty circle), |rfep from the adiabatic approximation (green filled diamond) and the numerical calculations with 
cut-off excitation number Cp = 2 (blue filled circle) and Cp — 3 (blue empty square). Here we consider the cases of ^ = 0.04a; 
and SI = cj, A = 0.04a; (a), Q, = 0.4a;, A = a; (b), SI = a;, A = 1.6a; (c) and SI = 0.4a;, A = 2a; (d). 



converge to fixed values and tk or the precise values 
of the single-photon scattering amplitudes. 



C. Results and discussions 

In Figs. 3 and 4, we illustrate the single-photon scat- 
tering amplitude rt given by our GRWA approach with 
the first and second order approximations, i.e., r'j}'' and 

(2) . 

rj. , the one given by the rotating-wave approximation 
and the result from the numerical diagonalization of the 
rotated Hamiltonian U^^HU. In our numerical calcula- 
tions the total excitation defined in p2p is cut off at a 
given number Cp for the Hamiltonian U~^HU, and the 
results with Cp = 2, 3 are shown in our figures. 

In Fig. 3 we calculate the refiection rate |rfcp. It is 

clearly shown that the results jj'fc^^P and I'r^.^^P from the 
first and second order approximation for (f/w) consist 
very well with each other. Therefore, in most of the cases 
with 1^1 << oj, the first order approximation is good 
enough for our GRWA approach. 

Furthermore, it is shown that in the case of Fig. 3(a) 
where the weak-coupling and near-resonance conditions 



are satisfied, both the results from the rotating-wave 
approximation and our GRWA approach fit well with 
the numerical calculations. Nevertheless, in the cases of 
Figs. 3(b)-3(d) where the rotating-wave approximation 
is not applicable, the results from our GRWA approach 
also consist significantly well with the numerical calcu- 
lations. This observation is further confirmed by Fig. 4 
where the real and imaginary parts of given by differ- 
ent approaches are illustrated. 

In Figs. 3(b) and 3(d) with ft = 0.4a;, we also com- 
pare our results with the ones given by the adiabatic 
approximation. It is shown that as we argued in Sec. 
II, the adiabatic approximation may be not applicable 
even when il << u, while our GRWA approach can still 
provide reasonable results. 

Therefore, the results in Figs. 3 and 4 show that, our 
GRWA approach with the first-order approximation for 
i^/uj), or our results r[.^^ and i^^"* in Eqs. ((36l [37| can be 
used as a good analytical approximation for the single- 
photon scattering amplitudes in the parameter region 
with 1^1 << a; and fl < uj. 



8 



n=e);X=S,=OM(a n=0.4(o;X=(o;^=0.04(o n=a);X=1 .6to;5=0.04(B £2=0.4cb;;i=2cb;^=0.04cb 




-0.3 -0.15 0.15 0.3 -0.3 -0.15 0.15 0.3 -0.3 -0.15 0.15 0.3 -0.3 -0.15 0.15 0.3 
k/(27t) k/(27i) k/(27t) k/(27t) 



n=(i);X.=^=0.04(fl n=0.4(o;X=(o;^=0.04(o n=(o;X=1 .6(o;^=0.04(o £l=0.4(o;X=2(o;^=0.04(o 




-0.3 -0.15 0.15 0.3 -0.3 -0.15 0.15 0.3 -0.3 -0.15 0.15 0.3 -0.3 -0.15 0.15 0.3 
k/(27t) k/(27i) k/(27t) k/(27i) 



FIG. 4: (Color online) The real parts (a-d) and imaginary parts (e-h) of the single-photon scattering amplitude rk given 
by our GRWA approach with the first-order approximation for (^/tij), i.e., r^^' in EQ. p6|) (red solid line), the rotating- wave 
approximation (green empty circle) and the numerical calculations with cut-off excitation number Cp = 2 (blue filled circle) 
and Cp = 3 (blue empty square). Here we consider the cases of ^ = 0.04a; and Q = uu, A = 0.04lj (a, e), Q — OAuu, X = uu {h, 
i),n = uj, A = 1.6uj (c, g) and Q = 0.4a;, A = 2cj (d, h). 



IV. THE SCATTERING AMPLITUDES IN THE 
STRONG-COUPLING CASE 

In the above section we derived the single-photon 
scattering amplitudes with our GRWA approach. It is 
pointed out that, our results in Eqs. p2all32el) are appli- 
cable for arbitrary large coupling between the TLS and 
the photon. Now we consider a special case where the 
TLS is strongly coupled to the photon in the resonator 
array, so that the condition 



A2 



« 1 



(38) 



is satisfied. We further assume the frequency il of the 
TLS is equal to or smaller than the photon frequency 
uj, i.e., < oj. In this strong-coupling case the expres- 
sions in Eqs. (|32all32el) can be significantly simplified and 
then one obtains simple pictures for both the quantitative 
calculation and the qualitative estimation of the single- 
photon scattering amplitudes. 

Under the condition ([38)) . we only keep the leading 

term proportional to (A^/a;^) exp (— 2A^/a;^) in uj, 



(0) 



and J(°'i) defined in Eqs. (I32aj|32ep . Then we have 



.(0) 

■'is 

(1) 
Iff 



2nx^ f A2 

— — cxp -2— 
1^ \ 



,(0) 



?,,(0) 

0. 



(39a) 

(39b) 
(39c) 



Therefore, in Ea. (|3T|) of the Hamiltonian H^^, we only 
need to keep the first two terms and the terms propor- 
tional to uj'fg and w^^'. This simplification implies that, 
in the strong-coupling regime our system is equivalent to 
the simple ID resonator array with the frequency of the 
0-th resonator shifted from w to w -I- , while the pho- 
tonic hopping intensity between the 0-th resonator and 
the ±l-st ones shifted from f to ^ -f Lo'^^g ■ Since we have 
also assumed < w, it is apparent that w^^^ << in 

the strong-coupling regime. Then the shift wj^'' of the 
photonic hopping intensity is negligible. We only need 
to consider the effect given by the frequency shift wj""* of 
the photon in the 0-th resonator. Namely, our system is 
finally equivalent to a ID resonator array, in which the 
0-th resonator has the frequency uj + Lo'fg , while all the 
other resonators have the same frequency uj. In this case 
the Hamiltonian H^^^^^ is approximated as 



H 



RWA(l) 
R 



4-, ^°^Jn 



(40) 



which leads to the single-photon scattering amplitudes 
rk 



"Iff 



J°J - 2i$sin(fc)' 



tk ^ rk + 1. 



(41) 
(42) 
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FIG. 5: (Color online) The single-photon reflection rate \rkf 
in the cases with strong TLS-photon coupling. Here we show 
the results given by Ea. (|41|l (red solid line), rotating- wave 
approximation (green empty circle) and the numerical cal- 
culations with cut-off excitation number Cp = 2 (blue filled 
circle) and Cp = 3 (blue empty square) for the cases with SI = 
CO, X = 1.6a;, ^ = 0.002cj (a), Q = 0.4a;, A = 1.4a;, ^ = 0.03a; 
(b), = oj, A = 2a;, ^ — 0.04a; (c). In the three cases we 
have a;jg'/5 — 15, 1, 0.06 respectively. 



ton would be reflected. Namely, in such a limit we have 
rfc « 1, tk ~ 0. Likewise, if the effective frequency shift 

Lo^iJ is much smaller than 4^, the frequency of the 0-th 
resonator would be approximately the same as the other 
resonators, and then every photon transmits through the 
0-th resonator. In this limit we have w and tfe w 1. 

In Fig. 5 we plot the photon reflection rate jrfep in the 
strong-coupling cases and compare the results given by 
our result in Eq. (j4ip and the one from numerical diago- 
nalization of the rotated Hamiltonian U~'^HU with cut- 
off excitation number Cp = 2, 3 respectively. It is clearly 
shown that our results in Eq. (|4T|) fit well with the nu- 
merical results. Furthermore, it is illustrated that in the 
case of Fig. 5(a) where we have H = 15, the photon 
reflection rate jr^p is almost unit for all the incident mo- 
mentum k. In the case of Fig. 5(c) with tofj /S, = 0.06, 
we have jr^ p « in the region with non-zero momentum 
k. All these observations are consistent with our above 
qualitative analysis. 

In the end of this section, we remark that, since all the 
quantities defined in Eqs. (|32all32ep exponentially decay 
to zero with (A/w), for any given values of f2, w and 
^, when the TLS-photon coupling intensity A is large 
enough, we can always neglect all these parameters and 
approximate the Hamiltonian j^^^^(i) g^g ^j-ee Hamil- 
tonian He for the array of identical resonators. There- 
fore, when the TLS-photon coupling is strong enough, 
the photonic scattering effect becomes negligible and we 
have Tfe = 0, tfe = 1 for the photon with any incident 
momentum k. 



V. THE GRWA FOR RESONATOR ARRAY 
WITH MULTI TLSS 

In the above sections, we generalized the GRWA to the 
system with ID resonator array coupled to a single TLS, 
and calculated the single-photon scattering amplitudes 
with the GRWA approach. In the end of this paper, we 
extend the GRWA to more general cases with m two-level 
systems coupled to the resonator array. For simplicity, 
here we assume each TLS is individually located in a 
resonator. Then the Hamiltonian of the total system is 
written as 



Hm — He + Ham + H 



IM 



(43) 



with the Hamiltonian He of the resonator defined in Eq. 
([1]), the Hamiltonian Ham of all the TLSs given by 



A straightforward result given by the above expres- 
sions of the scattering amplitudes is that, when the effec- 
tive frequency shift Lofg of 0-th resonator is much larger 
than the band width 4^ of the free Hamiltonian He of 
the array of resonators with the same frequency lu, the 
0-th resonator would be far-off detuned with the pho- 
ton with any incident momentum k, and then every pho- 



O ™ 



(44) 



13=1 



and the interaction Hamiltonian Him defined as 

m 

Him = A^(t('^^ {<^l(i3) + <^c(p)) ■ 

13=1 



(45) 
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Without loss of generality, here we assume the /3-th TLS 
is located in the c (/?)-th resonator. 

In such a general system, we straightforwardly develop 
the GRWA approach with the unitary transformation 
procedure in Sec. II and Sec. III. To this end, we first 
write the Hamiltonian Hm as 



— H 



Ml 



with Hm\ and Hm2 defined as 



H, 



M2 



Hmi — He + Him, 
Hm2 — Ham- 



(46) 
(47) 



Then we should find a unitary operator Um which can 
eliminate the linear terms of (a^, a}^) in Hmi and satisfies 



U mHmiUm 



4"-^-^ (a]aj+i + /i.e.) - Cm (48) 



single-photon scattering amplitudes under the conditions 
1^1 << uj and 57 < w. It is shown that in comparison with 
the rotating-wave approximation, our GRWA approach 
can give good results in a much broader parameter region. 
Especially, in the far-off resonance case with i7 << w, the 
adiabatic approximation is no longer applicable for our 
current system, while our approach still works well. We 
also generalize our GRWA approach to the ID resonator 
array coupled to multi-TLSs. 

In this paper, we assume the resonators in the ID array 
are single-mode ones. However, in the resonators used in 
the experiments, there usually exist more than one pho- 
ton modes. Then in the cases with strong TLS-photon 
coupling the multi-mode effect may be important. Like- 
wise, it may be also necessary to go beyond the two- level 
approximation and include the higher excited states of 
the artifical atoms in the strong-coupling cases. These ef- 
fects will be discussed in our coming work for the calcula- 
tion of the photon scattering in an multi-modc-resonators 
array or a multi-mode wave guide beyond the rotating- 
wave approximation. 



with Ca/ a constant c- number. The analytical calculation 
of Um is obtained in appendix D. 

With the operator Um, we apply the unitary trans- 
formation to the total Hamiltonian H, and make the 
rotating-wave approximation to the transformed Hamil- 
tonian U^HmUm- Finally we perform the inverse uni- 
tary transformation. Then the GRWA Hamiltonian for 
the resonator array with multi TLSs is obtained as: 



H 



M 



Hmg 

( oo 

Um \ ^ \PMn {UmHmUm) Pmu 



u 



M ' 



(49) 



which is a direct generalization of the result in Eq. ([20]). 
Here Pmu is the projection operator to the eigen-space 
of the total excitation operator 



El 

0=1 



e) 



a]aj 



(50) 



] = -oo 



with respect to eigen- value n. It is straightforwardly ob- 
served that, in the case of $7 << a;, such an approach also 
includes the intra-band transitions which are missed in 
the adiabatic approximation. On the other hand, under 
the weak-coupling condition and near-resonance condi- 
tion, this approach converges to the rotating-wave ap- 
proximation. 



VI. CONCLUSIONS 

In conclusion we generalize the GRWA to the hybrid 
system of a ID single-mode resonator array coupled to 
a single TLS, and obtain the analytical results for the 
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Appendix A: GRWA for the TLS-coupled 
single-mode bosonic field 

In this appendix, we introduce the GRWA approach 
proposed by Irish in Ref. [STj in the view of adiabatic 
approximation. To this end, we begin with the sim- 
ple Rabi Hamiltonian for the TLS-coupled single-mode 
bosonic field: 

-ffRabi = wfltfl + ^o-^ + A(Tj;(a'l' -I- a). (Al) 

Here a and are the annihilation and creation operators 
of the single-mode bosonic field with frequency uj respec- 
tively. f2 is the energy level spacing between the excited 
state |e) and the ground state \g) of the TLS, and A is 
the coupling intensity between the TLS and the bosonic 
field. The Pauli operators CTz and are defined in Sec. 
II. 

Since the Hamiltonian i?Rabi does not have simple in- 
variable subspaces, the exact diagonalization of i?Rabi is 
rather complicated [6l| . However, Jaynes and Gumming 
showed that [63 |. under the near-resonance condition 



\uj - n\ « \uj + 

and weak coupling condition 

|A| << UJ, Q, 



(A2) 
(A3) 
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the term |e) (5I a''' + h.c would be safely neglected. Then 
the Hamiltonian i?Rabi is approximated as the Hjc 
which is defined as 

Hjc ^i^a^a + X{\e){g\a + h.c). (A4) 

That is the so called rotating-wave approximation. Af- 
ter this approximation, the Hamiltonian iJRabi becomes 
invariable in the two-dimensional subspaces spanned by 
the states \g,n) and \e,n — 1) for n —1,2,..., as well as 
the one-dimensional subspace spanned by I5, 0) and then 
can be diagonalized easily. 

For the convenience of our discussions on the GRWA, 
here we introduce the projection operators P^2^^ defined 
as 

pin) ^ r \0g) {Og\ , n = 

R'^bi ]^\n,g) {n,g\ + \n- l,e) {n- l,e\, n>l ' 

(A5) 

Then the Jaynes-Cumming Hamiltonian Hjc in Eq. (jA4[) 
is re-written as 

00 

Hjc = T.PnliHn.,^PnL (A6) 

71=0 

Now we introduce the GRWA approach, which is devel- 
oped as an analytical approximate method to diagonal- 
ize the Hamiltonian H in Eq. (IA1[) in a broad parameter 
region where the rotating-wave approximation could be 
either applicable or not. The GRWA is closely related to 
both the rotating-wave approximation and the adiabatic 
approximation [GS'-'gs'I for the TLS-coupled single-mode 
bosonic field ;57t - r59i] which is used in the case of far-off 
resonance 

n«Lo. (A7) 

Therefore, before introducing the GRWA, we first intro- 
duce the adiabatic approximation in the system of a TLS 
and a single- mode bosonic field [stI [58| . In such an ap- 
proximation, the bosonic field is considered to be the 
fast- varying part of the total system and the TLS is con- 
sidered as the slowly- varying part. Then the Hamiltonian 
^^Rabi is rewritten as 

-ffRabi = ^^Rabil + -ffRabi2, (A8) 

where 

-ffRabii = wa'l'a 4- Xaxia^ + a) (A9) 

is the self-Hamiltonian of the fast-varying part together 
with the interaction between the fast-varying part and 
the slowly- varying one, and 

-ffRabi2 = ^cr, (AlO) 

is the self-Hamiltonian of the slowly- varying part. 

The Hamiltonian i/Rabii is easily diagonalized with the 
eigenstates: 



and the relevant eigen-energies: 

En± =ujin- X^/uj^). (A12) 

Here |±) are the eigenstates of ax with eigen-values ±1 
and \n±) are defined as 

\n±) — exp[=F A/w(a1' - a)] \n) . (A13) 

In the Rabi Hamiltonian the states \a, n) and \a' , n') are 
coupled by the term -ffRabi2- 

The spirit of the adiabatic approximation is described 
as follows [63l - l65j : under the far-off resonance condition 
Q. « uj, the motion of fast- varying part or the bosonic 
field adiabatically follows the slowly-varying part or the 
TLS, and can be frozen in the adiabatic branches with 
fixed quantum number n, or the two-dimensional sub- 
spaces spanned by |-|-, n) and | — , n) for n = 1,2, .... Then 
we neglect the iJFiabi2-induced transitions between the 
states \a,n) and \a' ,n') with n ^ n' . Then the eigen- 
states and eigen-energies of H are approximated as: 

l*±.n) = ^(|+,n)±h,n)), (A14) 

and 

E±^n - ±^{n-\n+) + L^{n - X^/lo''). (A15) 
respectively. 

Now we introduce the GRWA. In the "adiabatic basis" 
„)}, the Hamiltonian H is rewritten as 

i?Rabi = 5I (^Rabi):'„"'|*a,„)(*a',„'| (Al6) 

n,n' Q,a' — ± 

with the matrix elements 

(-ffRabi)"';,"' = («'a,n|-ffRabi|*a',„')- (Al7) 

In the GRWA, the Rabi Hamiltonian i?Rabi is approxi- 
mated as Hrq which is defined as 

-l-oo 

= EE(^R-bi):;:i*a,„)(*.,„i 

n— Q — di 

00 

+ 51 (^Rabi) + ;r' |*-,„)(*+,„-l| +/I.C.. 
n=l 

(A18) 

Namely, only the matrix elements (iJRabi)^ n inside 
each two-dimensional subspaces spanned by the states 
and |^'+^„_i) as well as the one in the one- 
dimensional subspace spanned by the state |^'-,o) are 
kept in the GRWA. In other words, in the GRWA one 
takes into account only the quantum transitions between 
the states |^'a,n) and l^* with 



|±,n) = |±)(g)|n±), 



(All) 



(A19) 
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with Na = 1 for a = + and Na = for a = — . 
Then, similarly as in the rotating-wave approximation, 
the Hamiltonian is reduced into a 2 x 2 blocked diagonal 
matrix in the GRWA. 

It can be shown that, under the weak-coupling and 
near-resonance conditions, the GRWA returns to the 
rotating-wave approximation. On the other hand, un- 
der the far-off resonance condition, the results given by 
the GRWA converges to the one from adiabatic approx- 
imation. Therefore, the GRWA smoothly connects the 
ordinary rotating wave approximation and the adiabatic 
approximation, and then works well in a more broad pa- 
rameter regime, especially the region with strong TLS- 
photon coupling and far-off resonant 51 < a; (see, e.g., 
Ref. 131 and Fig. 8 of Ref. [g^). 

In the end of this appendix, we emphasis that, the 
Hamiltonian Hug defined in (jA18l) could also be re- 
written as [52 



H 



RG 




Ur = exp 



P„ (iJ^^i/RabiC/fl) Pn) \ U^' (A20) 



(A21) 



with the unitary transformation Uji defined as 



Comparing Eq. (jA20l) and Eq. (jA6p , one can find that the 
GRWA is nothing but the "rotating-wave approximation 
for the rotated Hamiltonian U^^^HjxsihiUR" . 



Appendix B: The unitary transformation U 

In this appendix, we calculate the unitary transforma- 
tion U determined by Eq. (jH]) and Eq. ([TU)) . It is obvious 
that U can be expressed as the product of the displace- 
ment operators for each resonator mode, i.e., we have 



+ 00 



(Bl) 



J = -oo 

Then the Hamiltonian Hi is transformed into 



U-^HiU 



+ [ujao - £,{ai + + A] a^ial + gq) - C. 

(B2) 

The exact expression of C is not required here. Compar- 
ing the result Eq. (|B2I) with Eq. we get the equations 
for the parameters {aj}: 

Lja, -CK+i = 0,j =±1,±2,... (B3) 
wao - + a-i) + A = 0,j = 0. (B4) 



Now we solve Eqs. (jB31 IB4p with two steps. First, we 
introduce a cut-off for the equations at j = ±A^, and 
solve the equations 



AN) 



(B5) 

^(4 



(N) 



A = 0, J = 0, 



a 



{N) 
N 



a 



(N) 

-N 



(B6) 
(B7) 



A straightforward calculation shows that 



JN) 



UJl 



AN) 



2C. 



with UJl given by 



(B8) 



(B9) 



Therefore, under the weak- hopping assumption |f | << w, 
we have ^/wi < 1. Substituting Eq. (|B7)) into Eq. (|B8t . 
we obtain 



(AT) 

"O = 



-Xuji 



{^f\N\ _ 1 



and then we obtain all from Eq. (jBSp . 

Second, we consider the solutions of Eqs. (IB5IIB7P in 
limit of TV — >■ oo as a trial solution of Eqs. (|B31 IB4[) : 



N\ 



(BIO) 



tto = lim a, 



(N) 



Xloi 



2e 



OJUtl OJl 



(Bll) 



Substituting Eq. (|BTT1) into Eqs. dBHEH, we find that 
the latter ones are exactly satisfied. Therefore, {aj} from 
Eq. (|B11|) are the realistic solutions of Eqs. ([Ill EH. 
Then the unitary transformation U defined in Eq. ([8]) 
takes the form Eq. ^ with aj given by Eq. ([TU|) . 

We emphasis that, as shown in Eq. (|Bip . U is the prod- 
uct of the displacement operators for each resonator, the 
magnitude of the displacement for the mode in the j- 
th resonator is described by \aj\. Furthermore, since 
^/wi < 1, the result Eq. (IB11[) implies that \aj\ expo- 
nentially decays with Then we have 



lim exp [aj da; (a] — aj)] = 1 



(B12) 



Therefore, for the resonators which are far away from the 
TLS, the relevant displacements in U would be negligible, 
this is consistent with our above considerations. 
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Appendix C: The proof of Eq. (fT9l) 



In this appendix, we prove Eq. (jl9[) in our maintext. 
To this end, we first define the operator Uk [Pk] and the 
state \n (fc)) for the ID resonator array as a function of 
the number Pk- 

Uk m = exp {-2/3fc [a (k)^ - A (fc)] } (CI) 



respectively. Here |0) is the vacuum state of the resoator 
array, A (k)^ is defined in Eq. ([6]). We further define the 
function /(/3) as 



/(/3) = n('^W \l^kWk]\n (k)) 



with 



and 



\nik)) = ^l=Aik f'^'^ |0) 



(C2) 



fW) = Y[{n{k)\Uk[f3k]\n'{k)) 

k 

= []e-|2/^-l'/2(n(fc)|e-2/^'=^W'e2/3,A«|n'(fc)) 



n 



n'(fe) 

E 

m=Max[0,n'(A:)— n(fc)] 



n (fc) |yl (fc)^™ A (fc)" |n' (/c)) 



ml [n (fc) — n' (fc) + m]! 



(C3) 



/3 = (/3fe,,/3fe,,...). (C4) 
The straightforward calculation shows that 



(C5) 



which gives 



/(-/3) = (-l)^(")-^("') /(;3) 



(C6) 



for N{n) = N {n') . (C7) 



On the other hand, with the above definitions and 
straightforward calculations, we have 

^ /(/3o); (C8) 

- /(-/3o), (C9) 

where |±,n) are the eigen-states of the Hamiltonian Hi 

defined in Sec. IIB and /3o — (/3ofei , /3ofe2 7 •■•■)• Here we 
have 



Pak — 



1 



'2tt 



+ 00 

E 



aie 



ikl 



(CIO) 



with ai defined in Eq. ([T0|. Then using Eq. ([T4| and 
Eq. (fT7|). wc re- write the matrix element H^;'^_ as 



(* + ,ff|i^2|*-,«') 

^{{+,n\-,n')-{-,n\+,n')) 

n 

1 



fWo) - fhf3o) 



(Cll) 



Therefore, our above result in Eq. (IC7p directly leads to 
Eq. (HH). 



Appendix D: The unitary operator Um 

In this appendix, we calculate the unitary operator 
Um defined in Similar as in appendix B, it is easy 

to prove that Utm is the product of the displacement 
operators of each resonator mode: 



+00 



(Dl) 



To derive the expression of aj , we define the column 
vector 






(D2) 



Then the straightforward calculation shows that, the con- 
dition (pS]) is equivalent to the equation 



ujd — S^Ka — A. 



(D3) 



Here K is, a square matrix with the element Kij in the 
i-th row and ^'-th column given by 



(D4) 



In Eq. (|D3|) . A is a constant column vector with the j-th 
component Aj defined as 



/9=1 



(D5) 
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Therefore, we formally have the expression of a: 



(D6) 



Furthermore, we notice that the matrix w — is di- 
agonalized as 

uj-^K^ J dk{uj-2^ cos k) v{k)v{k)^ (D7) 

with the j-th component of the column vector v(k) sat- 
isfies 



Then we have 



1 



dk 



v{k)v{k)^ 



(D8) 



(D9) 



oj — S^K J w — 2^ cos k 
Substituting (|D9I) into (|D6p , we get the expression of a*^ 



a 



M _ ^ ^ / dfc£^^Pl!%_£M. (Dio) 



2tt 



/3=1 



w — 2^ cos k 



In the case of a single TLS located the 0-th resonator, we 
have m = 1 and c (1) = 0. Then af^ is expressed as 



exp [ifcj] 
w — 2^ cos k 



(Dll) 



On the other hand, in such a single-TLS case, the value 
of a*^ is also given by (|Blip . Therefore, we have 



A f dh ^^p['^^'] - 

27r 7 w — 2^cosfc 2^^ — woji 



(^)l^l (D12) 



with wi defined in Sec. II. B. Substituting (ID12[) into 
(jPlOp . we finally obtain 



, i \\i~-cm 



-(-)'■ 



(D13) 



Therefore, we get the analytical expression of the unitary 
operator Um defined in Eqs. (pSllDlj) . 
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